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Abstract. The dimensional reduction of D- dimensional spacetimes aris- 
ing in string/M-theory, to the conformal Einstein frame, may give rise to 
cosmologies with accelerated expansion. Through a complete analysis of 
the dynamics of doubly warped product spacetimes, in terms of scale in- 
variant variables, it is demonstrated that for D > 10, eternally accelerating 
4-dimensional k = — 1 Friedmann cosmologies arise from dimensional reduc- 
tion on an internal space with negative Einstein geometry. 



1. Introduction 

The current standard model of cosmology has as an essential element the ac- 
celerated expansion of the universe. In order to achieve accelerated expansion 
the strong energy condition must be violated. Numerous matter models which 
provide accelerated expansion have been proposed, among them the cosmolog- 
ical constant A and scalar field models like quintessence and /j-essence. 

It has been shown by Townsend and Wohlfarth [10] by considering a Kaluza- 
Klein reduction of a D-dimensional spacetime with a hyperbolic internal space 
that the dimensionally reduced universe may exhibit a period of accelerated 
expansion, even though the Z}-dimensional model one starts from does not 
violate the strong energy condition. This result thus circumvents a 'no-go' 
theorem which states that dimensional reduction of supergravity models arising 
from string/M-theory (which in particular satisfy the strong energy condition) 
cannot give rise to cosmologies with accelerated expansion in case the internal 
geometry is time-independent [7, 9]. 

The D-dimensional models discussed in [10] are vacuum spacetimes which 
are warped products of a (3 + l)-dimensional k = Friedmann model with an 
n-dimensional hyperbolic space, and are thus examples of Lorentzian doubly 
warped products. The dimensionally reduced model exhibit a transient phase 
of acceleration. In [1] it has been shown that it is possible to obtain late time 
accelerated expansion, if the (3 + l)-dimensional model is a k = — 1 Friedmann 
model. In the same paper it was also argued, based on a perturbation argument, 
that solutions with eternal acceleration might exist. 

In this paper we apply the method of scale invariant dynamics to the study of 
vacuum doubly warped products and their dimensional reduction. In particular, 
we prove the existence of eternally accelerating solutions for D > 10. Our 
analysis shows that eternally accelerating cosmologies cannot be obtained from 
perturbations of the (dimensionally reduced) D-dimensional Friedmann model, 
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as was argued in [1]. For D < 10 we prove that no eternally accelerating models 
can occur. 

We consider line elements onaD = m|n|l dimensional spacetime M. x 
M x N of the form 

-dt 2 + a 2 (t)g + b 2 (t)h , (1.1) 

where (M, g), (N, h), are m- and n-dimensional Einstein spaces of non-negative 
Einstein curvature. (Positive curvature is discussed in Sec. 6.) The vacuum 
Einstein equations give a system of ODEs for a, b; however, in our approach, 
we formulate the equations in scale invariant variables and obtain a regular 
dynamical system on a compact state space, cf. Sec. 2. The state space is topo- 
logically a disk; in the interior there are represented models where both factors 
g, h have negative Einstein curvature; solutions on the boundaries correspond 
to models where one of the factors is Ricci flat. 

For the D-dimensional geometry, we show that the behavior is asymptoti- 
cally Friedmann in the expanding direction, corresponding to the approach to 
a stable fixed point (F*) in the formulation of the scale invariant dynamics, 
and asymptotically Kasner-like in the collapsing direction. The asymptotic be- 
havior in the expanding direction depends on the dimension; for D < 10, the 
Friedmann point (F*) is a stable spiral point, while for D > 10, (F*) is a stable 
node. The system exhibits two qualitatively different types of orbits for the 
collapsing direction: the generic orbit is asymptotically Kasner like with non- 
vanishing generalized Kasner exponents, which implies that in the direction of 
the singularity one of the factors in the doubly warped product is expanding 
while the other is collapsing. This case is analogous to the vacuum Kasner 
spacetime with line element —dt 2 + t 2p dx 2 + t 2g dy 2 + t 2r dz 2 and with exponents 
(p,q,r) = (2/3,2/3,-1/3). In addition, there are two exceptional orbits, as- 
ymptotic to one of two fixed points (Fa), (Fb) on the boundary of the state 
space, with the property that one of the factors converges to a constant and 
the other scales asymptotically with proper time toward the singularity. These 
cases are analogous to the flat Kasner spacetime with exponents (1,0,0). The 
analysis leading to this description is performed in Section 3, in particular, 
Theorem 3.1 summarizes the main results. 

The issue of whether the D-dimensional models considered in this paper are 
stable in the sense that they admit large families of perturbations with quiescent 
behavior at the singularity is discussed in Section 4. We find that for D > 11, 
this is indeed the case. This agrees with the general analysis of [6, 3]. Thus we 
see an indication in this simple context that spacetimes with D = 11 exhibits 
special features. 

In Section 5 we apply the results about the L>-dimensional scale invariant 
dynamics to the dynamics of the dimensionally reduced models. We take the 
factor (TV, h) as internal space. It turns out that a cosmological model that 
arises through dimensional reduction exhibits accelerated expansion when the 
corresponding solution of the scale invariant dynamical system lies in a certain 
region A of the state space, cf. Fig. 3. This domain of acceleration A inter- 
sects the boundary of the state space in an open interval; hence, solutions on 
the boundary (which describe models where g is Ricci flat) exhibit a phase of 
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acceleration, which explains the transient acceleration found by Townsend and 
Wohlfarth [10]. Furthermore we find that the Friedmann fixed point (F*) lies 
on the boundary of A. The fact that, for D < 10, the Friedmann point (F*) is a 
spiral point thus entails that models cannot exhibit late time accelerated expan- 
sion, even with negatively curved spatial slices and hyperbolic internal space. 
Instead one finds that cosmologies which arise through dimensional reduction 
of spacetimes with D < 10 have an infinite sequence of episodes with acceler- 
ating and decelerating expansion. However, for D > 10, based on arguments 
combining the local properties of (F*) and global properties of scale invariant 
dynamical system, we are able prove that there exists a unique dimensionally 
reduced cosmology with eternal acceleration (i.e., a cosmology that expands at a 
accelerated rate for all times), cf. Theorem 5.1. Note, however, that uniqueness 
of this model amounts to the model being unstable under any perturbation; the 
behavior of generic models is described in Corollary 5.2. 

Finally, in Section 6 we consider the case where one of the factor in (1.1) has 
positive Einstein curvature and one negative curvature. Using again a formu- 
lation in terms of scale invariant variables, we are able to perform a complete 
analysis of the dynamics of the D-dimensional geometries as well as their di- 
mensional reduction. The domain of acceleration can be found explicitly, and it 
turns out that there exist models with one or two epochs of expansion (whereof 
one can be partly accelerating), but that all models undergo recollapse to a big 
crunch. 

2. Preliminaries 

Let (M, g) and (N, h) be Einstein manifolds of dimension m and n with 

Ric 9 = k g (m + n — l)g , Ric^ = kh{m + n — l)h ; (2.1) 

kg and k^ are constants taking values in {+1, 0, —1}. On the D = 1 + m + n 
dimensional spacetime R x M x N consider a line element of the form of a 
doubly warped product 

-dt 2 + a 2 (t)g + b 2 (t)h, (2.2) 

where a > and b > without loss of generality, and impose the vacuum 
Einstein equations. Let i,j be indices running over 1 . . . m + n. The second 
fundamental form is 

R 3 = | -(«/«) <V 1 ^ hJ ^ m 

\ — (b/b) 6/ m + l<i,j<m + n, 

where we use the notation / = dtf. Let p = —a/a, q = —b/b. Then the mean 
curvature is H = trK = mp + nq, and the Einstein evolution equations imply 

p = pH + k g m + n ~ 1 (2.3a) 

or 

, fn + n — 1 , , 

q = qH + k h - 2 . (2.3b) 

From these evolution equations we obtain 

H = mp 2 + nq 2 . 
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Here we have employed the Hamiltonian constraint equation, 

, m(m + n — 1) , n(m + n — 1) „ 9 . 9 9 . 

= k g ^ = + k h — -= + H — (mp 2 + nq 2 ) . 

a 2 b l 

Note that k g < and kh < entails / 0; in the following we focus on this 

case, the case k g > 0, kh < is treated in Sec. 6. We now introduce scale 

invariant variables (P,Q, A, B) according to 

The variables A, B are curvature quantities. By definition we obtain mP+nQ = 
1, and the Hamiltonian constraint now reads 

1 = -{jn + n - l)(/e 9 m,4 2 + k h nB 2 ) + (mP 2 + nQ 2 ) . 

In the following we introduce the time t by d T = H~ l dt, we will use a prime ' 
to denote differentiation w.r.t. r. Note that with our conventions H < 0, hence 
by introducing the time r we have the singularity to the future. 



3. Analysis 

The variable transformation (a, b,p, q) ^ (A, B, P, Q, H) enables us to write 
the Einstein equations as a system of evolution equations, 

H' = H(mP 2 + nQ 2 ) (3.1) 

and 

A' = A [P - (mP 2 + nQ 2 )} (3.2a) 

B' = B[Q - (mP 2 + nQ 2 )} (3.2b) 

P' = P [1 - ( m p 2 + n Q 2 )} + (m + n- l)k g A 2 (3.2c) 

Q' = Q [1 - ( m p 2 + nQ 2 )} + (m + n- l)k h B 2 , (3.2d) 

supplemented by two constraint equations, 

d = mP + nQ - 1 = (3.3a) 

C 2 = (mP 2 + nQ 2 ) -(m + n- l)(k g mA 2 + fc^riB 2 ) -1 = 0. (3.3b) 

Since Eq. (3.1) decouples, the entire dynamics is encoded in the reduced dynam- 
ical system (3.2), which is autonomous and regular for all (A, B,P,Q) 6 R 4 . 
From (3.2a) and (3.2b) it follows that A > and B > are invariant under the 
flow of the system; henceforth, without loss of generality, we will always impose 
these conditions, i.e., we consider the dynamical system (3.2) on the state space 
X = {(A, B, P, Q) | (A > 0) A (B > 0)}. 

In our subsequent analysis it will turn out that there exist solutions with 
AB —>■ for \t\ —>■ oo. This suggests to include the boundaries A = and 
B = of the state space in our dynamical systems analysis; note that the 
system (3.2) can be smoothly extended to X. 

The equations on the invariant subset A = (respectively B = 0) can be 
interpreted as the reduced system of coupled equations that arises when the 
first factor (respectively the second factor) of the metric (2.2) is Ricci flat. 
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Figure 1. The state space S. 



This is because setting A = in (3.2) and (3.3) corresponds to setting k g = 
and discarding the decoupled equation for A. (In the case k g = 0, due to the 
decoupling, the equation for A does not carry any dynamical information; note, 
however, that the equation must be added in order to reconstruct the original 
variables.) 

By construction the constraints (3.3) are preserved during the evolution. To 
see the propagation of constraints explicitly we compute 

C[ = [l- (mP 2 + nQ 2 )} d - C 2 (3.4a) 
C 2 = -2(mP 2 + nQ 2 ) C 2 . (3.4b) 

Hence the physical state space 

S := {(A,B,P,Q) | (Ci = 0) A (C 2 = 0) A (A > 0) A (B > 0)} (3.5) 

and the closure S are invariant subsets of X (and X, respectively). 

When the variable constraint (3.3a) is solved for Q, the Hamiltonian con- 
straint (3.3b) becomes 

— (m + n)(p — Y -(m + n-l)(k Q mA 2 + k h nB 2 ) = m + n ~ . ( 3 6 ) 

n V m + nJ a m + n 

When kg = kh = —1, this condition defines an ellipsoid centered at (A,B,P) = 
(0,0, l/(m + n)). Hence, S corresponds to a quarter ellipsoid, which is topo- 
logically a disk; see Fig. 1. In particular, we see that in the case kg — k^ — 1 
the state space S has a compact closure. (Recall that the dynamical system 
induced on the boundaries of S represents the reduced dynamics of the cases 
kg = 0, kh = —1 and k g = —1, k^ = 0.) 

In the subsequent sections, k g = = — 1 is understood (the case kgkji — 
— 1 is treated in Sec. 6). This entails that the dimensions of the factors M and 
N satisfy m > 1 and n > 1. 
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3.1. Equilibrium points and invariant subsets. An elementary analysis 
shows that the equilibrium points in S of the autonomous system (3.2) are 

(Fi )2 ) : A = B = 0, mP 2 + nQ 2 = 1, mP + nQ = 1 , (3.7a) 

(F A ) : A = 0, B 2 = n ~ 1 1 -, P = 0, Q = -, (3.7b) 

n z m + n — 1 n 

(F B ) : A 2 = m ~ 2 1 1 -, 5 = 0, P=-, Q = 0, (3.7c) 

to 2 m + n — 1 m 

(F») : A = B = P = Q = — *— . (3.7d) 

m + n 

Note that in (3.7a) the equations mP 2 + nQ 2 = 1 and mP + nQ = 1 possess 
two solutions (Pi,Qi) and (P 2 ,Q 2 ), which define the fixed points Fi and F2, 
respectively; we have Pi > (Q\ < 0) and P 2 < (Q 2 > 0); see Fig. 2. Note 
that \P{\ < 1 and \Qi\ < 1, because m > 1 and n > 1. 

As noted above, the boundaries of the state space are invariant subsets. We 
denote U A = {A = 0} n S and U B = {B = 0} n S. The set U A contains the 
fixed points Fi, F 2 , and Fa; Ub contains Fi, F 2 , and F^. Recall that these 
invariant subsets describe the dynamics of a spacetime with spatial geometry 
which is a product of an Einstein manifold with a Ricci flat space. 



3.2. Stability of the fixed points. 

The case (Fi i2 ): Since the fixed points (Fi ;2 ) on S are given as intersections 
of the one-dimensional invariant subspaces Ua and Ub, the eigenvectors of 
the linearization of the system (3.2) must be tangential to Ua and Ub- For 
(Fi), the associated eigenvalues are (Pj — 1) and (Qi — 1), which follows from 
j 4 _1 A'|(f.) = Pi — 1 and the analogous relation for B~ 1 B'. Since |Pj| < 1 and 
\Qi\ < 1, the eigenvalues are negative, and we conclude that the points (Fi i2 ) 
are sinks. 

The case (F a ,b): As the fixed point (Fa) lies on Ua, one eigenvector of 
the linearization of the system at (Fa) is tangential to Ua- The associated 
eigenvalue is (n — l)/n, which follows when we set A = in (3.2c) and compute 
P _1 P'|(F A ) = 1 — (1/ro). There exists a second eigenvector which is transversal 
to Ua; the associated eigenvalue is (— 1/n), since j 4~ 1 j 4'|(f a ) = — (1/n). We 
conclude that (Fa) is a saddle, so that generic orbits in S do not tend to (Fa)- 
However, there exists exactly one orbit that converges to (Fa) along the stable 
subspace as r — > 00. Finally, since Ua coincides with the unstable manifold, 
(Fa) is a repellor in Ua- A similar analysis applies to (Fb). 
The case (F*): Let J* denote the linearization matrix of the system (3.2) 
at (F#), i.e., the Jacobian of the right hand side of (3.2). The tangent space 
T( F „)S of the constraint manifold S at (F*) in X is spanned by the vectors 
(n,-m,0,0)*,(0,0,n,-m)*. Let 



S = 



( 


n 


^ 




-m 










n 


V 





—m J 



ETERNAL ACCELERATION FROM M-THEORY 



7 



and denote by S 1 any (2x4)-matrix that is left-inverse w.r.t. S. The restriction 
of J* to T(f,)S is given by S~~ 1 J*S, which yields 

^*It ( f*)S ^ m + n (— 2(m + n — 1) m + n — l) ' 

The eigenvalues of this matrix are 

y = m + n-1 ± y^m + n- l)(m + w-9) 

We distinguish three qualitatively different cases: when (m+n) < 9, there exists 
a non- vanishing imaginary part; when (m + n) = 9, X\ = X^; when (m + n) > 9, 
the eigenvalues are real and 1 > A^ > A£ > 0. For all (m,n), the real part of 
A^ 2 is positive, which entails that (F*) is a repellor, i.e., locally stable toward 
the past. 

3.3. Global dynamics. The global properties of the flow on the boundaries of 
S are simple: since Ua and Ub are one-dimensional, the local stability analysis 
implies the global dynamics. Ua can be viewed as an interval, whose end points 
are the stable fixed points (Fi) and (F2). In between there exist one additional 
fixed point, the repellor (Fa). It follows that every orbit (different from the 
fixed points) originates from (Fa) and ends in (Fi) or (F2). The picture on Ub 
is analogous. 

On S, we observe that P is monotone when P < 0, i.e., P' < when P < 0, 
which is because mP 2 + nQ 2 < 1. Analogously, Q' < when Q < (which in 
turn corresponds to P' > when P > (1/m) via (3.3a)). The flow on S is thus 
particularly simple for P < and P > (1/m). 

Now consider the non-negative function Z = A m B n on X (Z is essentially 
the inverse of the rescaled volume density of the spatial metric). Under the 
side-condition C2 = 0, cf. (3.3b), hence in particular on S, Z attains a global 
maximum, namely at the point (F*). On the boundaries Ua and Ub, Z be- 
comes zero, which is its minimal value. From (3.2) we obtain that 

Z' = Z[l- (m + n)(mP 2 + nQ 2 )] (3.10) 

on S. Since (mP 2 + nQ 2 ) > l/(m + n) on S unless P = Q = l/(m + n), the 
bracket in (3.10) is negative almost everywhere. In the special case P = Q = 
I /(m + n) we obtain Z" = and 



Z'" = -2Z(m + n)(m + n-l) 



>n(A 2 -(^) 2 ) 2 +n(B 2 -(^) 2 ) 



(3-11) 

which is non-positive and vanishes only when A = B = l/(m + n), i.e., at the 
fixed point (F*). We infer that Z is strictly monotonically decreasing along all 
orbits in S\(F*). This excludes that there exist any nontrivial periodic orbits 
in S and allows us to invoke the monotonicity principle: the a-limit of every 
orbit in S must coincide with the fixed point (F*); the w-limit of every orbit 
must lie on the union of Ua and Ub, i.e., AB — > as r — > 00. 
Summarizing, we have proved the following theorem: 
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(Fi) 



(Fi) 



(Fb) 



w 




(Fa) 



(Fa) 



(F 2 ) 



(a) m + n < 9 (D < 10) 



(b) m + n > 9 (D > 10) 



Figure 2. Schematic of the flow on the state space S for the 
cases m + n < 9 and m + n > 9. The point at the center is the 
fixed point (F*). (The flow in the case m + n = 9 looks qualita- 
tively like (b); however, the eigenvectors of the linearization at 
(F*) coincide.) 

Theorem 3.1 (Global dynamics). Let k g = kh = —1 andm,n > 1 and consider 
the dynamical system (3.2) on the two-dimensional physical state space S. 

(1) Consider an orbit with A = 0, B / that is different from (Fa). The 
a-limit of the orbit is the fixed point (Fa), the oj-limit is (Fi) or (F2). 
The statement for A 7^ ; B = is analogous. 

(2) Consider the family of orbits with AB 7^ that are different from (F*). 
The a-limit of every orbit is the fixed point (F*). There is one orbit 
whose uj-limit is (Fa), another one whose uj-limit is (Fb). The lj -limit 
of a generic orbit is one of the equilibrium points (Fi), (F2). 

The flow of the dynamical system on S is depicted in Fig. 2. 

The behavior of the decoupled variable H is given by Eq. (3.1): the equation 
implies that H is monotonic in r. By virtue of the inequalities 1 > (mP 2 + 
nQ 2 ) > l/(m+n) on S we obtain H < H' < H/ (m+n); in particular, H \ —00 
asr/ +00, and H / as r \ —00. 

Remark. In [4, 5] the D-dimensional Einstein vacuum equations for (Riemann- 
ian) metrics dt 2 + a 2 g + b 2 h are analyzed as a constrained Hamiltonian system, 
in particular as regards the integrability of the equations. The authors provide 
evidence that the cases D = 10, 11 are special (and perhaps integrable); indeed, 
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at least for three subcases conserved quantities for the flow are constructed 
in [4], see also the remark in Sec. 5. 

3.4. Asymptotics of the scale factors a and b. 

Every orbit in S gives rise to a solution —dt 2 + a 2 g + b 2 h of the D-dimensional 
Einstein vacuum equations via the transformation (2.4). The particular asymp- 
totic behavior of the orbit translates to characteristic asymptotic behavior of 
the scale factors a and b: 

Proposition 3.2 (Asymptotic behavior). An orbit that converges to (Fi^) as 
t — > oo generates a solution with asymptotic behavior of Kasner type, i.e., 
a ~ t p , b ~ t q as t \ 0, where mp + nq = 1 and mp 2 + nq 2 = 1. When an orbit 
converges to (F*) ost-> — oo, then the corresponding solution is of Friedmann 
type as t — > oo: a ~ t, 6 ~ t. Finally, convergence to (Fa,b) as t ^ oo leads to 
a — > const, 6 ~ i and a ~ t, 6 — > const as i \ 0, respectively. 

Remark. Note that the orbits on the boundaries converge to (Fa,b) as t ^ — oo. 
The corresponding solutions satisfy a — ► const, 6 ~ i and a ~ i, 6 — > const as 
t — > oo. 

In the following we establish the above proposition; furthermore we analyze 
in more detail the approach to Friedmann as t — ► oo. 

The case (Fi^): Consider a solution that converges to (Fi) as r — ► oo. We 
have -ff — > — oo, which suggests that the limit corresponds to a singularity. 
Since dt/dr = H^ 1 and H decreases exponentially with r, we find that r — > oo 
corresponds to i \ 0. 

In a neighborhood of the fixed point (Fi) the dynamical system (3.2) on S 
can be approximated by its linearization at (Fi), which is given by 

A! = A[P 1 -1], B , = B[Q 1 -1], P = Pi, Q = Qi, (3.12) 

cf. the previous stability analysis. It is important to note that the equa- 
tions (3.12) coincide with the system of equations (3.2) & (3.3) obtained for 
kg = kh = 0. We conclude that all solutions that converge to (Fi) as r — ► oo 
behave asymptotically like solutions representing a spacetime where both fac- 
tors g and h are Ricci flat. In other words, the dynamical effects of curvature 
become negligible at the singularity. 

When we solve (3.12) and the decoupled equation (3.1) for H, and recall that 
dt/dr = H -1 , we obtain via (2.4): 

a = a t Pl , b = b t Ql (3.13) 

for t \ 0; ao and &o are constants. By construction, mP\ + nQ\ = 1 and 
mP 2 + nQ 2 = 1; Pi > 0. A completely analogous result is obtained for solutions 
that approach (F2), where (Pi,Qi) is replaced by {P2,Q2}- 
The case (F*): Consider a solution that converges to (F*) as r — > —00. We 
have H —>■ and t —>■ 00, which suggests that (F*) corresponds to the limit of 
infinite expansion. Since P = Q in the limit we expect Friedmann like behavior. 

Assume that m + n < 9. In a neighborhood of the fixed point (F*) the 
dynamical system (3.2) on S is approximated by its linearization (3.8). Using 
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the variable transformation 
B 



1 y/m + n — 1 



y/9 - (m + n) V° 4 V™ + n - 1 
in T( Ft )S, the linearized system takes a normal form given by 

yjm + n 



(3.14) 



B 



_1 ( J *It (f , )S ) 



B 



yjm + n - 
2(m + n 



1/ \/m + n — 1 
,-\/ 9 ~ (m + n) 



(m + n 



V9 

\An + n — 1 



Based on this we can show that, modulo terms of order 0(t ( m+n 2 )), 



(3.15) 



a = t 



1 - nc y/9 -{m + n) t ^( m+n ' 1 )/ 2 cos(-) + 

+ nc y/m + n - lt-^ rn+n - 1 ^ 2 sin(-) 
1 - mc y/m + n-lt- (m+n - 1 ^ 2 sin(-) + 



(3.16a) 



(3.16b) 



+ mc y/9 -(m + n) t-( m +™- 1 )/ 2 

cos • , 

where the argument of the trigonometric functions is 

(•) = (-±Vm + n-\yj9-(m + n) log(i/d )) ; (3.17) 

Co, do are constants. 

When m + n > 10 the eigenvalues AJ 2 of the linearization of the dynamical 
system at (F*) are real with 1 > AJ > 2A 2 > 0. Oscillatory terms as in (3.16) 
do not occur in this case. A thorough analysis of the system yields that orbits 
converging to (F*) as r — > — oo correspond to solutions of the type 



a 



t 

b = t 



l + nc t- x " +0(t- 2X ") 



1 - mc t' x * + 0(t 



,-2A* 



(3.18a) 
(3.18b) 



as t — > oo, where A* = (m + njX^ = \(m + n — 1 — y/Jm + n — l)(m + n — 9)). 
(In (3.18) there appears only one free parameter; the second parameter of the 
two-parameter family of solutions is connected to terms of the order £-( m + n ) A i .) 

In the special case m + n = 9 the linearization of the dynamical system at 
(F*) is represented by a non-diagonalizable matrix. Orbits converging to (F*) 
as r — ► — oo correspond to solutions of the type 

a = t [1 + nco(logt)^ 4 + nd t~ A + 0((logt) 2 t" 8 )] (3.19a) 

b = t [1 - mc (logt)t~ 4 - md t~ 4 + O ({log t) 2 t~ 8 )] (3.19b) 

as t — > oo, where Co, do are constants. 

The case (Fa,b) : Consider a solution along the orbit that converges to (Fa) 
as r — > oo. We have H — > — oo and t — > as r — > oo, which suggests that the 
limit corresponds to a singularity. Using again approximation techniques we 
eventually find that 

b= Vm jHZ 1 t (3.20) 



as t 



yjm + n — 1 

ao b = — t 

y/n — 1 

0. This completes the proof of Proposition 3.2. 
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Corollary. Combining the above results with the statements on the behavior of 
generic orbits in S, cf. Theorem 3.1, we find that generic solutions of the D- 
dimensional vacuum equations are of Kasner type, as in (3.13), as t — ► and of 
Friedmann type, as in (3.16), (3.18), or (3.19), as t — ► oo. Interestingly enough, 
the approach to the Friedmann solution is oscillatory if and only if m + n < 9 
{D < 10). 

4. AVTD CONDITION 

In this section we briefly address the question of whether there exist families 
of perturbations (that are general in the sense that they depend on the maximal 
number of free functions) of the .D-dimensional models described above that 
exhibit quiescent behavior at the singularity, cf. [6, 3]. 

Consider a (d + l)-dimensional Kasner spacetime, with Kasner exponents pi, 
i = 1 . . . d, and line element 

d 

-dt 2 +J2 t2pi ( dxi f ■ 

i=l 

The vacuum Einstein equations imply the Kasner relations 

i=l i=l 

We consider the case where d = m + n, with pi = P, i = 1, . . . m, and pi = Q, 
i = m + 1 . . . m + n. Then the Kasner relations (4.2) read 

mP + nQ = 1, mP 2 + nQ 2 = 1 , (4.3) 

which we recognize as the constraint equations (3.3) in case A = B = 0. These 
equations hold at the equilibrium points (Fi^), cf. (3.7a). 

Following [6], see also [3, §3], the condition for asymptotically velocity-term 
dominated (AVTD) behavior for an asymptotically Kasner spacetime is 

1 + Pi ~ Pd ~ Pd-i > , 

which in terms of the generalized Kasner exponents P, Q reads 

1 + P - 2Q > . (4.4) 

At the equilibrium points (Fi^), we have 



1 / m + n-l\ 1 / jm + n 

Pi,2 = ■ — 1 ± nd , Qi >2 = ■ — 1 =F my 



m + n \ V ran } ' m + n \ V mn 

1.5) 

The condition (4.4) at (Fi), i.e., 1 + Pi — 2Qi > 0, is satisfied for all (m,n) 
such that 

An 

m > , n > 1 . (4.6) 

n — 1 

Equivalently, 1 + P 2 — 2Q 2 > 0, when m > 1 and n > 4m/(m - 1). Eq. (4.6) 
implies, with d = m + n, 

d > ^±^> („ > 1) . (4.7) 
n — 1 
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For n = 2, 3, 4, 5, . . . the right hand side is 10, 9, 9|, 10, . . ., and for n > 3, the 
right hand side is monotonically increasing. Therefore we find, in agreement 
with the result of [6], that d > 10 is a necessary condition for AVTD behavior. 



5. Accelerating cosmologies from compactification 

Dimensional (Kaluza-Klein) reduction transforms classes of D-dimensional 
vacuum spacetimes into classes of (1 + m)-dimensional spacetimes (where typ- 
ically m = 3) with nonlinear scalar fields. In some instances the resulting 
(1 + m)-dimensional models represent accelerating cosmologies, as has been 
demonstrated by a sizable number of examples, see, e.g., [10, 11, 12] and [1, 2]. 
In the following we present a systematic treatment of accelerating cosmologies 
arising from the compactification of the D-dimensional spacetimes discussed 
in Sec. 3. In particular, we prove existence and uniqueness of a model that 
exhibits eternal acceleration. 

Consider the family of D = 1 + m + n dimensional vacuum spacetimes 
M x M x N with metrics — dt 2 + a 2 g + b 2 h that has been constructed. We 
perform a dimensional reduction, i.e., we assume that the m spatial dimensions 
connected to the factor g are the large spatial dimensions representing the clas- 
sical spacetime (where typically m = 3 in order to obtain a four-dimensional 
spacetime), while the n spatial dimensions connected to h are to be compacti- 
fied; 

-fl? + a 2 g + b 2 h reduction ^ = ^^2 + Q 2 g) ({u) 

We choose a conformal Einstein frame for the 1+m metric, see, e.g., [2, App. A], 
i.e., we conformally rescale —dt 2 + a 2 g by choosing <p according to 

e 4> = b -n/(m-l) j (5.2) 

which entails that the reduction of the (l + m + n)-dimensional Einstein-Hilbert 
action results in an effective action, again of Einstein-Hilbert type, for the metric 
7 and a nonlinear scalar field ip that stems from the compactified dimensions, 
tp = -(87r) _1 /2 n -i/2(- m _ i)i/2 ( m + n _ !)i/2 0_ Therefore, (M,^,(p) satisfies 

the Einstein nonlinear scalar field equations, i.e., 

Ric 7 — iScal 7 7 = 8irT , (5-3) 
where T is the energy- momentum tensor of the scalar field, 



liiv (//, v = . . . m) . (5.4) 



The potential V(ip) of the scalar field is an exponential function, 



t 7- / \ , Q x-i n{m + n- 1) 
1/(99) = (»7r) — 1 exp 



.m+n—1 

8ir-==\ p 

yjm - 1 V n 



(5.5) 



For details on dimensional (Kaluza-Klein) reduction see, e.g., [8]. 

For the metric 7, by introducing a new time variable i through di/dt = 
^n/(m-i) we f] na i}y obtain 

7 = 6 2n /(™-l)(_^2 + a 2 g ^ = _ d p + ^n/(«-l) X2 ^ = _ dt "2 + ^ 
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with a = ab n /( m ~ l \ 

In the following we investigate the spacetimes R x M with metrics 7; in 
particular we focus on the question of whether these solutions give rise to cos- 
mological models that exhibit accelerated expansion. 

By construction, every orbit in S is associated with a solution 7 = — dt 2 + 
d 2 g of the (1 + m)-dimensional Einstein equations (with nonlinear scalar field 
<p), where the scale factor a = afe™^" 1-1 ) is determined via (2.4). A simple 
calculation shows that the derivative of the scale factor a is 

§ = ^>-i- <"> 

Since P < 1 on the entire state space S, da/dt is positive for all solutions, so 
that all metrics 7 describe expanding cosmologies. 
Differentiating (5.7) we obtain 

d 2 n nh _n /( m_1 ) r 1 

^4 = H 2 -(l-P) 2 -(m+n-l) (kJm- l)A 2 + k h nB 2 ) . (5.8) 

dt z m — 1 L 

The first term in brackets is negative, the second term is non-negative on the 
state space S. Clearly, in a neighborhood of the fixed points (Fi^), where 
A = B = 0, the sum is negative, hence d 2 ajdt 2 < 0. However, in the following 
we will establish the existence of a domain ACS such that 

d 2 a 

> for all (A, B, P, Q) G A C S . (5.9) 

dt z 

The domain A is the domain of acceleration: whenever an orbit passes through 
this domain, the cosmological model (R x M, 7) it represents undergoes accel- 
erated expansion. 

Remark. As before, in (5.7), k g = kh = — 1 is understood. However, recall that 
setting kg = 0/kh = corresponds to setting ^4 = 0/B = 0; thus the cases 
kg = 0, kh = — 1 and k g = — 1, kh = are included in our treatment by our 
investigation of the boundaries of S. 

Consider first the boundary Ua given by A = in S. Setting A = in (5.8) 
we find by simple algebraic manipulations of the bracket in (5.8) and by using 
the Hamiltonian constraint that 

d 2 d m + n-l k mn ( m + n\ 2 

n~V~~k — (5 - 10) 

with k mn = m(m + n) + n. We infer that d 2 d/dt 2 > in a neighborhood of 
P = (m + n)/k mn , namely for P G (7r_, 7r+), where 



m + n ± ny/m + n — 1 
vr± = ; (5.11) 

Kmn 

for all m and n, ir± are positive. Hence, on the boundary A = there exists 
a non-empty interval of acceleration; note that the fixed point (Fa) is not an 
element of this interval. 

Remark. Since the orbit (Fa)-(Fi) on the boundary Ua passes through an 
interval of acceleration, the corresponding model (R x M, 7) exhibits a (finite) 
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(Fi) (Fi) 




(F 2 ) (F 2 ) 
(a) m,n arbitrary (b) m + n > 9 (D > 10) 

Figure 3. The domain of acceleration A on the state space S. 
Whenever an orbit passes through this domain, the cosmological 
model (M x M, 7) it represents undergoes accelerated expansion. 
For m + n > 9 (m > 2) there exists a unique orbit that lies 
entirely in A; the associated model exhibits eternal acceleration. 



phase of accelerated expansion; note that M is Ricci flat for this model. This 
solution was originally found in [10]; it provided the first example for accelerated 
expansion from compactification. 

On the boundary Ub, by setting B = in (5.8) and using the constraint, we 
obtain 

7P«-> + »-»( p -k)' ■ (5 - 12) 

Hence, cPa/dt 2 = at the fixed point (Fb) and negative elsewhere on B = 0. 

Inserting the Hamiltonian constraint into Eq. (5.8) it is straightforward to 
show that dA (given by cPa/di 2 = 0) defines an ellipse in S. Combining this fact 
with the results collected above we conclude that the domain of acceleration A 
in S possesses the following main properties: A intersects the boundary ^4 = 
in an interval that does not contain (Fa) and the boundary B = in the point 
(Fb); dA is tangential to the boundary B = at that point. Moreover, a simple 
calculation shows that (F*) lies on dA; see Fig. 3(a). 

The boundary of A consists of three semipermeable segments: the flow of 
the dynamical system on the segment (^4 = 0) — (F*) is directed toward the 
interior of A, while the flow on the segments (F*) — (Fb) and (Fb) — (^4 = 0) 
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is directed toward the exterior. This is because 
d fd 2 a\ ( „ 1 \ 2 



8t \ dt 2 



cx - P [P . (5.13) 

<Pa/dP=o \ m J \ m + n; 



Consider the case m + n < 9 (D < 10). The cosmological model (R x M, 7) 
possesses a phase of accelerated expansion whenever the associated orbit in 
S passes through A. Superimposing Fig. 2(a) and Fig. 3(a) we see that a 
generic solution (represented by an orbit converging to (Fi) or (F2)) does not 
exhibit accelerated expansion for early times. However, toward the future, there 
exist infinitely many phases of accelerated expansion (associated with the orbit 
spiraling out from (F*)). The properties of the non-generic solution converging 
to (Fa) are analogous, but the solution converging to (Fb) is of a different 
kind: it generates the single model that is accelerating for early times (i.e., in 
a neighborhood of the singularity). 

Now consider the case m + n > 9 (D > 10), i.e., Fig. 2(b) together with 
Fig. 3(a). The figures suggest that, in contrast to the case m + n < 9, there 
might exist one orbit that lies entirely in A. We thus formulate the following 

Theorem 5.1 (Existence and uniqueness of eternal acceleration) . Letm+n > 9 
(with (m,n) 7^ (2, 7)); k g = kh = —1. Then there exists a unique solution 7 = 
—dP + a 2 g of the (1 + m)- dimensional Einstein equations (with nonlinear scalar 
field (p) arising from the dimensional reduction of a (1 + m + n)- dimensional 
vacuum solution, such that d 2 a/dt 2 > for all t. 

Proof. Since (Fi), (F2), and (Fa) do not lie in A, orbits converging to any of 
these fixed points cannot generate models with eternal acceleration. This leaves 
the orbit connecting (F*) with (Fb) as the only candidate. To establish eternal 
acceleration for the associated model we must show that this orbit lies entirely 
in A. To that end consider the intersection of the plane 



(A,B,P,Q) eX cosV' M -sin^ [P =0^ , (5.14) 

\ m + nj \ m + nj J 

where ip is a constant, and the state space S C X. This intersection generates 
a curve a$ in S that connects (F*) with a point on B = 0, cf. Fig. 3(b). When 
ip is chosen according to 

ip = arccot (^y/m + n—l [y/m + n— 1 + y/m + n — 9]^ , (5.15) 

then ay, possesses the following favorable properties: (i) connects (F*) with a 
point on B = with P > 1/m, i.e., with a point on the (Fb) — (Fi) segment (and 
does not intersect the segment (F*) — (Fb) of dA); (ii) is semipermeable 
for the flow of the dynamical system: orbits can pass through in one direction 
only (namely from the lower left side to the upper right); see Fig. 3(b). We have 
thus constructed a region A^, C A, enclosed by and segments of dA, that is 
past invariant for the flow of the dynamical system. Since the orbit converging 
to (Fb) lies in A^ as r — > 00 it follows that it must lie in A^ and thus in A for 
all times. By this the theorem is established. 
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Remark. For one special case of m and n, namely m = 2, n = 7, the above argu- 
ment fails (in particular (ii) does not hold) . A thorough numerical investigation 
yields that this is indeed because the statement is wrong: the orbit connecting 
(F*) with (Fb) does not lie entirely in A but intersects <9A transversally on the 
segment (A = 0) — (F*); hence the orbit lies in A for early times and outside 
for late times. The numerical evidence thus suggests that in the special case 
(m, n) = (2, 7) eternal acceleration is impossible. 

Remark. In three special cases the "orbit of eternal acceleration" can be ob- 
tained explicitly: for (to, n) G {(2, 8), (3, 6), (5, 5)} the orbit is given simply 
as 

1 1 \ / . 1 \ / . 1 \ /„ 1 



A — ) - ( A(f b) — I I P — I = 



K m m + n J \ m + n J \ v ' m + n J \ m + n y 

(5.16) 

where Ar F \ = (m — l) 1 / 2 m~ 1 (m + n — 1) _1//2 , cf. (3.7c). This is analogous 
to the analysis of [4], in which it is shown that the D-dimensional Einstein 
vacuum equations (for Riemannian metrics) regarded as a Hamiltonian system 
are integrable for (m, n) £ {(2, 8), (3, 6), (5, 5)}. This coincidence suggests that 
the analysis of [4, 5] might also apply in the Lorentzian case. 

Corollary 5.2. The generic solutions cannot exhibit eternal acceleration. Or- 
bits that converge to (Fi) correspond to solutions that exhibit a phase of ac- 
celerated expansion {t\,t ); there are two subcases: solutions corresponding to 
orbits closer to the boundary A = possess a finite phase of accelerated expan- 
sion, i.e., < t[ < t < oo, while solutions corresponding to orbits closer to 
B = accelerate forever toward the future, i.e., < t[ < t Q = oo. In contrast, 
accelerated expansion does not occur for all orbits that converge to (F2): orbits 
closer to A = generate solutions with decelerating expansion. (Note that the 
same holds for the non-generic solution converging to (Fa).^ The expansion for 
orbits closer to B = is accelerating in (t\,t ) with < t\ < t Q = 00. 

The proof of the corollary is based on the established properties of the orbit 
of eternal acceleration (F*)-(Fb); (m, n) g" {(2, 7), (7, 2)} is necessary. 

We conclude by discussing the leading term of the scale factor a(t) of the 
model (R x M, 7) as i — > and i — ► 00. 

(Fi 5 2): For the generic orbit that approaches the Kasner fixed point (Fi) 
or (F2) as t — > we obtain t oc t m ^ 1 ~ p ^^ n ~ 1 \ hence t —> corresponds 
to t — ► and we obtain a oc i 1 / m as i — > 0. 

(F*): In the limit we have t oc ^( m +"-- 1 )/( m - 1 ) an( j the leading term in a 
is proportional to i. 

(Fa,b): The orbit converging to (Fa) generates a solution that satisfies 
t oc t(m+n-i)/(m-i) an d a oc fV(™+«-i) i n the limit t 0. The solution 
corresponding to the orbit converging to (Fb) satisfies i oc t and a oc t 
in the limit t — ► 0. 

Remark. The fixed point (Fa) represents a model — dt 2 + a 2 g where g is Ricci 
flat; the scale factor satisfies a power law, a oc t n /( m + n ~ 1 ) . This model is 
well-known, it is the simplest of solutions of the Einstein nonlinear scalar field 
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equations with exponential potential (5.5). Since the exponent in V((p) is over- 
critical, the power in the scale factor a is less than one. 



6. The case k g > 0, k h < 

Consider again the metric —dt 2 + a 2 g + b 2 h with Ric g = k g (m + n — l)g 
and Ric/j = kh(m + n — l)h. Assume k g = 1 and kh = — 1. For the analysis 
of this case the equation system (3.2) is ill-adapted since the variables become 
ill-defined (which is due to the possibility of H going through zero). We replace 
the system (3.2) by equations that are adapted to the case k g = 1, kh = —1. 
With 

D=( H 2 +k g m{m+ j i - i) y (6.i) 

define, in close analogy to the previous definitions, 

P = j>, Q = ^, A=(aD)-\ B=(bD)- 1 . (6.2) 

This transforms the Einstein equations to a system consisting of the decoupled 
equation 

D' = D ((mP + nQ)(mP 2 + nQ 2 ) + k g m(m + n- l)PA 2 ) (6.3) 
and the coupled equations 

A' = A [P - (mP + nQ)(mP 2 + nQ 2 ) - k g VA 2 } (6.4a) 

B' = B[Q- (mP + nQ)(mP 2 + nQ 2 ) - k g VA 2 ] (6.4b) 

P' = P [{mP + nQ)(l - mP 2 - nQ 2 ) - k g VA 2 ] + k g (m + n - 1)A 2 (6.4c) 

Q' = Q [{jnP + nQ)(l - mP 2 - nQ 2 ) - k g VA 2 ] + k h {m + n - 1)B 2 , (6.4d) 

where k g VA 2 = k g m(m + n — 1)PA 2 and the prime denotes d T = D~ 1 dt- The 
two constraints are 

Ci = (mP + nQ) 2 + k g m(m + n- 1)A 2 -1 = (6.5a) 
C 2 = mP 2 + nQ 2 - k h n(m + n- 1)B 2 -1 = 0. (6.5b) 

We consider the system (6.4) on the (new) state space S defined by C\ = 0, 
C 2 = and A > 0, B > 0. In Eqs. (6.4) and (6.5) k g = 1 and k h = -1 is 
understood. Note that the dynamical system induced on the boundaries A = 
and B = of the state space describes the reduced dynamics of the cases k g = 0, 
kh = — 1 and k g = 1, k^ = 0, respectively. 

Equation (6.5b) together with B > describes the "northern hemisphere" 
of an ellipsoid. By (6.5a), A > implies — 1 < (mP + nQ) < 1. Hence, S 
corresponds to a strip on the ellipsoid that connects two diametrically opposite 
segments of the "equator" B = via the "north pole" P = Q = 0. Topo- 
logically, S is a rectangle, whose sides are B = (with P > 0), A = (with 
mP + nQ = 1), B = (with P < 0), and A = (with mP + nQ = -1); see 
Fig. 4. 
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6.1. Equilibrium points and invariant subsets. The equilibrium points in 
S of the autonomous system (6.4) are 

(F1...4) : A = B = 0, mP 2 + nQ 2 = 1, mP + nQ = ±1 , (6.6a) 

(F A+ ,_): A = 0, B 2 = ^—± -, P = 0, Q = ±-. (6.6b) 

n A m + n — 1 n 

The equations mP 2 + nQ 2 = 1 and mP + nQ = ±1 possess four solutions 
(Pi, Qi), ■ ■ ■ , (P4, Q4), which define the fixed points Fi, . . . , F4. For i = 1,2 we 
have mPi + nQi = 1, Pi > (Qi < 0) and P 2 < (Q 2 > 0); for i = 3, 4 we have 
mPj + nQ, = -1, see Fig. 4. Note that (P 3 ,<5 3 ) = -(Pl,Qi) and (Pa,Qa) = 
— (P2,Q2)- Clearly, |Pj| < 1 and \Qi\ < 1, because m > 1, n > 1. We define the 
boundaries of the state space Ua = {^4 = 0} n S and Ub = {B = 0} fl S. 

6.2. Stability analysis. 

The case (F1...4): Since the fixed points (F1...4) are given as intersections of 
the one-dimensional invariant subspaces Ua and Ub, the eigenvectors of the 
linearization of the system (3.2) must be tangential to Ua and Ub- The relation 
J 4 _1 A'|( Fi ) = Pi — sign(mPj + nQi) and the analogous relation for B gives the 
eigenvalues of the linearization of the dynamical system at the fixed points. 
The eigenvalues are (Pj — 1) and (Qi — 1) for i = 1,2; hence (Fi^) are sinks. 
The eigenvalues are (Pj + 1) and (Qi + 1) for i = 3,4; hence (F^^) are sources. 
The case (Fa + _ ) : One eigenvector of the linearization of the system at (Fa + _ ) 
is tangential to Ua- The associated eigenvalue is ±(n — l)/n, which follows 
when we set A = in (6.4c) and compute P" 1 P / | (Fa+ j = ±[1 - (1/n)]. The 

second eigenvector is transversal to Ua; the associated eigenvalue is (=Fl/n), 
since j 4 _1 t1 / |(-f a+ ) = =F(l/n). We conclude that (Fa + _) are saddles, so that 
generic orbits in S do not tend to (Fa + _). However, there exists exactly one 
orbit that converges to (Fa + ) along its stable subspace as r — ► 00, and there 
exists exactly one orbit that converges to (Fa_ ) along its unstable subspace as 
r — ► —00. 

6.3. Global dynamics. The existence of periodic orbits in S is excluded by 
the absence of interior fixed points; furthermore, there exist no heteroclinic 
cycles. Consequently, all orbits converge to the fixed points as r — > ±00. 

The fixed points (Fa + _) are of particular interest. Recall that there exists 
exactly one orbit that originates from (Fa_) into the interior of S, and one 
orbit that converges to (F A+ ) as r — > 00. Making use of global properties of 
the flow we prove that these orbits do not coincide: 

Consider Eq. (6.4c) in the case |P| <C \Q\, A < B, i.e., 

P » k g (m + n - l)A 2 + (-k h )n 2 (m + n- 1)B 2 PQ . (6.7) 

In a neighborhood of (Fa_), we obtain P' > if P < 0. It follows that 
the orbit originating from (Fa_) satisfies P(r) > for all r sufficiently small. 
Analogously, in a neighborhood of (Fa + ), P' > if P > 0, hence the orbit 
converging to (Fa + ) satisfies P(r) < for sufficiently large values of r. Since 
P' > when P = by (6.4c), it follows that the two orbits are distinct. 
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Q 

(F 2 ) 




m,P + nQ = — 1 



Figure 4. A projection of the state space S in the case k g = 1, 
kh = — 1 and schematic of the (projection of the) flow. 



Using again the fact that P' > when P = it is straightforward to con- 
clude that the orbit that originates from (Fa_) converges to (Fi) as r — ► oo; 
conversely there is one orbit that originates from (F3) and ends in (Fa + ). 

6.4. Asymptotics of the scale factors. 

The case (F1...4): Consider a solution that converges to the fixed point (F3) 
as r — > —00. In a neighborhood of the fixed point the dynamical system (6.4) 
can be approximated by its linearization at (F3), which is given by 

A' = A[P 3 + 1] , B' = B[Q 3 + 1], P = P 3 , Q = Qs- (6.8) 

The equations (6.8) coincide with the system of equations (6.4) & (6.5) obtained 
for kg = kh = 0. Consequently, all solutions that converge to (F3) as r -> 
— 00 behave asymptotically like solutions representing a spacetime where both 
factors g and h are Ricci flat. In other words, the dynamical effects of curvature 
become negligible in the limit. 

Solving (6.8) and the decoupled equation (6.3) we obtain 

a = a t Pl , b = b t Ql (6.9) 

for t \ 0; ao and 60 are constants. We have used here that — P3 = P\ and 
—Q3 = Qi', recall that mP\ + nQ\ = 1, mP^ + nQ\ = 1, Pi > 0. A completely 
analogous result is obtained for solutions that approach (F4), where (P\,Qi) is 
replaced by (P 2 , C^)- 

Similarly, solutions that tend to (Fi) asr^oo correspond to 

a = a (T-t) p \ b = b (T - t) Ql ; (6.10) 

the analogous result holds for (F2). In (6.10), T is the time of the future 
singularity of the spacetime. 
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The case (Fa + _): Now consider a solution along the orbit that converges to 
(Fa + ) as r — > oo. We find 

a^a b= " /m t^ 1 (T-t) (6.11) 
yn — 1 

ast^T. Analogously, 



Vm+n-1 

a — ► ao o = ■== — t (6-12) 

yn — 1 

for a solution along the orbit that converges to (Fa_) as r — > — oo. 

The global properties of the flow of the dynamical system translate to the 
following statements: all solutions undergo recollapse; there exist two generic 
types of past singularities and one special type; analogously, there exist two 
generic types of future singularities and one special type, cf. Fig. 4. (I) Every 
solution with a past singularity of the type (a oc t Pz , b oc t^ 2 ) possesses a future 
singularity of the type (a oc (T — t ) Pl , b oc (T — t)® 1 ). (II) For a solution with 
a past singularity of the type (a oc t Pl , b oc t® 1 ) there exist three scenarios: 
(1) (a oc (T - t) p \ b oc (T - t)Qi), (2) (a oc (T — t) p \ b oc (T - t)Q 2 ), or (3) 
(a — > const, b oc (T — t)) as i — > T. The scenario (3) is not generic. (Ill) There 
exists a non-generic solution with a past singularity of the type (a — > const, 
b oc i) for t — ► 0. This solution possesses a future singularity of the type 
(a oc (T-t) Pl , 6oc (T-t)^ 1 ). 

6.5. Accelerating cosmologies from compactification. 

Consider, in complete analogy to Section 5, the dimensionally reduced space- 
time R x M with metric 

1 = -dt 2 + a 2 g, (6.13) 

where dt/dt = 6 n /( m_1 ) and a = ab n ^ m ~ 1 \ Orbits in S generate cosmological 
models (R x M, 7); in the following we analyze their properties. 

A simple calculation shows that the derivative of the scale factor a is 

hence the equation [(m — 1)P + nQ] = divides the state space: there exists a 
domain of expansion and a domain of contraction in S; the former is the region 
below the straight dashed line in Fig. 5. As long as an orbit lies in this domain, 
the cosmological model it represents is expanding. Superimposing Figs. 4 and 5 
we see that all models (R x M, 7) are expanding initially and contracting toward 
the end. 

Differentiating (6.14) we obtain 



d 2 a ab- n /( m -V „ r 1 
-D z 



-P 2 - — (m + n)Q 2 ] . (6.15) 
m m J 



dt 2 m — 1 

It is straightforward to show that d 2 a/dt 2 is positive in some domain A of the 
state space S. This domain is indicated by curved dashed lines in Fig. 5. The 
intersection with the domain of expansion is non-empty, thus we see that there 
exists a domain Ae C S of accelerated expansion: whenever a solution passes 
through this domain, the cosmological model (R x M, 7) it represents undergoes 
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Figure 5. The state space S in the case k g = 1, kh = — 1. 



accelerated expansion. (In the complementary domain Ac we have da/dt < 
and d 2 a/dP > 0.) 

The boundary of Ae consists of semipermeable segments. This is because 
^r(^£] oc -(1 - mP 2 - nQ 2 )(P - Q) (6.16) 

and 

A(*f\ oc(m-l)A 2 -ni? 2 . (6.17) 

dr\dt J da/dt=o K ; y ' 

Hence, orbits can enter Ae through the upper or the right boundary and leave 
through the lower boundary. 

We conclude that there exists a one-parameter family of solutions that ex- 
hibits a phase of accelerated expansion, cf. Figs. 4 and 5; these solutions are a 
subfamily of the family of solutions that connects (F3) with (Fi). Note that 
there can be at most one phase of accelerated expansion. Interestingly enough, 
not all solutions are simply expanding-contracting; there exist solutions with 
two phases of expansion, i.e., models that are expanding-contracting-expanding- 
contracting. For those models the second expanding phase begins as a phase of 
accelerated expansion. 



7. Concluding remarks 

In this paper we have analyzed the dynamics of D-dimensional vacuum dou- 
bly warped product spacetimes and their dimensional reduction to the confor- 
mal Einstein frame. The method of scale invariant dynamics has allowed us 
to give a comprehensive description of the global dynamics; in particular, we 
have proved that, for D > 10, when the 4-dimensional spacetime is a k = — 1 
Friedmann model and the internal space a negatively curved Einstein space, 
there exists a unique (thus non-generic) dimensionally reduced model which is 
expanding at an acclerated rate for all times (eternal acceleration). Generic 
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models, on the other hand, fall into classes: there exist models that exhibit a 
transient phase of acceleration, late time acceleration, or no accelerated expan- 
sion at all. 

The D-dimensional models we have considered are of the simplest type: dou- 
bly warped product solutions of the Einstein vacuum equations. Generalizations 
include multiply warped product spacetimes, see, e.g., [1], D-dimensional grav- 
ity coupled to n-form field strengths [15] or n-form field strengths and a dilaton, 
see, e.g., [13, 14]; for reviews see [11, 12]. 

In the study of these models, in particular in view of the question of accel- 
erated expansion, dynamical systems methods and phase space analysis have 
proved to be powerful tools; for instance, the methods permit the complete clas- 
sification of the qualitative late-time bahavior of models [11, 15] (as re-derived 
by the local dynamical systems analysis in Section 3.4 and at the end of Sec- 
tion 5). In the present article, in order to obtain a description of the global 
dynamics of models we have employed the method of scale-invariant dynamics. 
This approach naturally embeds the local analysis (late time behavior, . . . ) into 
a global context, which has allowed us, in particular, to establish the existence 
of a non-generic model exhibiting eternal acceleration. Whether the methods 
generalize to models as in [1, 13, 15] is under current investigation. 
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